Measurement of the polar-angle distribution of leptons from 
W boson decay as a function of the W transverse 
momentum in pp collisions at a/s=1-8 TeV 
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Abstract 

We present a measurement of the coefficient 0:2 of the leptonic polar-angle distribution 
from W boson decays, as a function of the W transverse momentum. The measurement 
uses an 80ib4 pb~^ sample of pp collisions at ■\/s=1.8 TeV collected by the CDF detector 
and includes data from both the W ^ e + u and W ^ fi + u decay channels. We fit the W 
boson transverse mass distribution to a set of templates from a Monte Carlo event generator 
and detector simulation in several ranges of the W transverse momentum. The measure- 
ment agrees with the Standard Model expectation, whereby the ratio of longitudinally to 
transversely polarized W bosons, in the Collins-Soper W rest frame, increases with the W 
transverse momentum at a rate of approximately 15% per 10 GeV/c. 



PACS numbers: 12.15.Ji, 12.38.Qk, 13.85.Qk, 13.38.Be, 14.70.Fm 



I. INTRODUCTION 



According to the Standard Model (SM), the polarization of W bosons produced at high 
transverse momentum (p^) is strongly affected by initial-state gluon radiation and quark- 
gluon scattering (the QCD leading-order diagrams for high-pr W production are shown in 
Fig. 1). The angular distribution of the leptons from the W —>■ £ + u decay reflects the 
changes in the W polarization. In the Collins-Soper W rest frame [1] the dependence of the 
cross section on the leptonic polar-angle at hadron level can be parametrized as 

-r— — oc (1 - Qai cos Ocs + "2 cos^ Ocs), (1) 
dcos6'c5 

where Q is the lepton charge. The effects of QCD contribute to the coefficients ai and 02, 
which are functions ofpj^ . Fig. 2 shows the theoretical expectation for ai{pjf) and a2{pj^), 
neglecting a correction from sea-quarks, calculated up to next-to-leading order in QCD [3, 4]. 
Sea quarks give an opposite sign contribution to the cos 6cs term when the W is produced 
by an antiquark in the proton and a quark in the antiproton, reducing the value of ai . Only 
in the limit pj^ — GeV/c, when ai = 2 and a2 = 1, does Eq. (1) describe the distribution 
of leptons from a transversely polarized W boson: da/dcos6cs oc (1 — Q cos Ocs)^, which 
is typical of a pure V — A interaction. As a2 decreases, the contribution from longitudinally 
polarized W bosons increases and so does the probability for the decay lepton to be emitted 
at large polar angle. On the other hand, ai measures the forward-backward leptonic-decay 
asymmetry. Fig. 2 indicates that the asymmetry is reduced at higher p^ . 

Understanding how QCD corrections affect lepton angular distributions is important 
in the measurement of the W mass {Mw) and rapidity distributions in pp experiments. 
The lepton angular distribution changes the shape of the transverse mass distribution, 
which is used to measure Mw- It has been estimated that an uncertainty of ±1% on 0:2 
corresponds to a shift of the measured value of M^y, determined by fitting the transverse 
mass distribution, of approximately ±10 MeV/c^ [5]. 

We present the measurement of a2 at various W transverse momenta, using both the 
electron and muon channels. The sensitivity for a measurement of ai is too low, due 
to the fact that the sign of cos Ocs is undetermined without a full reconstruction of the 
kinematics of the neutrino from the W decay. Hence, the only sensitivity to a\ comes from 
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FIG. 1: The QCD leading-order processes that give rise to W production at high-p^. In the 
top diagrams a giuon is radiated from one of the scattering quarks. In the bottom diagrams a 
quark-gluon scattering produces a W, together with a quark. 



the correlation between the geometrical acceptance of the detector and the phase space of 
the observed events. The current best measurement of 02 is reported in [4]. The results 
presented here reduce the uncertainty on 02 by about 50% up to pj!^ ~ 30 GeV/c, and are 
of comparable uncertainty at higher transverse momenta of the W. 

For completeness, the cross section differential in the azimuthal and polar lepton angles 
can be expressed in the most general form as 

dV 3 dV™^ 



d{pj^y dy dcosOcs d4>cs 16vr d{p}pydy 



[(1 + cos" 9cs) + 
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+^^0(1 - 3 cos^ 0cs) - QM sin 10 cs cos 0cs + 
+^^2 sin^ dcs cos 24>cs + ^3 sin Qcs cos ^cs + 
-QA4 cos dcs + M sin^ ^cs sin 24)cs + 
-Q sin 2Qcs sin 4>cs + M sin dcs sin 4)cs\ , 



(2) 



where y is the rapidity of the W boson, a'^'^'^ is the total (angle integrated) rate, and the 



2 



Ai terms weight the relative contributions to the total cross section due to the different 
polarizations of the W boson. By integrating Eq. (2) over (f) and comparing with Eq. (1) it 
follows that 

2^4 2 - 3^0 

"^ = 2Ti^' "^ = ^t:47' 

which relates the ai and a2 with the coefficients. The A^ coefficients are explicitly 
calculated in [3, 6]. 

This paper is structured as follows: Sections II and III describe the CDF detector and 
the W boson data sample, Sections IV and V outline the measurement method and detail 
the Monte Carlo event generator and detector simulation. Section VI contains the estimate 
of the background to the W data sample and Section VII summarizes the fits and the sys- 
tematic uncertainties. The results and conclusions are presented in Section VIII. 



II. THE COLLIDER DETECTOR AT FERMILAB (CDF) 

A complete description of the CDF detector can be found elsewhere [7] . We describe here 
only the components relevant to this work. CDF uses a cylindrical coordinate system (r, <p, 
z) with the origin at the center of the detector and the z-axis along the nominal direction 
of the proton beam. We define the polar angle 9 as the angle measured with respect to 

the z-axis and the pseudo-rapidity (77) as 77 = — ln(tan(^/2)). A schematic drawing of one 
quadrant of the CDF detector is shown in Fig. 3. 



A: Tracking 

The CDF tracking system in Run I consists of three tracking detectors: a silicon vertex 
detector (SVX'), a vertex time projection chamber (VTX) and an open-cell multiwire drift 
chamber (CTC). The tracking system is immersed in a 1.4 T solenoidal magnetic field 
aligned with the 2;-axis. The silicon vertex detector [8] is a silicon microstrip detector that 
covers a region in radius from 2.86 to 7.87 cm. It is divided into two identical "barrels" 
which surround the beampipe on opposite sides of the z = plane. Each barrel consists 
of four radial layers of silicon strip detectors, and each layer is divided in azimuth into 30° 
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FIG. 2: Theoretical NLO-QCD calculation of and a\ vs. p^. The limit GeV/c is the 

Quark Parton Model, for which q;2 = 1 and ai = 2 . 

wedges. The microstrips run parallel to the z direction so that the SVX' tracks particles 
in r — (/). The VTX [9] is a set of 28 time projection chambers, each 9.4 cm in length, 
surrounding the SVX' detector. It provides the z position of the interaction point with 
a resolution of 1 to 2 mm. The CTC [10], which extends out to a radius of 138 cm and 
|z| < 160 cm, measures a three-dimensional track by providing up to 60 axial and 24 stereo 
position measurements. The basic drift cell has a line of 12 sense wires strung parallel to 
the z-axis for axial measurements or 6 sense wires tilted ±3° in for stereo measurements. 
The set of all drift cells located at the same radius from the origin of the detector is called 
a super-layer. 

In this analysis the CTC is used for the tracking and VTX and SVX' are only used to 
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provide vertex information. The CTC track is constrained to point to the event vertex. 
The z location of the vertex is determined with the VTX, and the position in r — ^ is 
determined from the beam hne measured with the SVX'. The result of this procedure is a 
significant improvement in the CTC resolution. The momentum resolution of such tracks 
is u{pt)/pt = [(0.0009 prf + (0.0066)2] V2 ^-^^-^ pr measured in units of GeV/c. 

B: Calorimetry 

The CDF calorimetry is provided by four different calorimeter systems with a nearly con- 
tiguous coverage out to \r]\ = 4.2 . Three of the four systems have both electromagnetic (EM) 
and hadronic (HA) calorimetry. They are called "Central" (CEM, CHA), "Wall" (WHA), 
"Plug" (PEM, PHA) and "Forward" (FEM, FHA). The central and wall calorimeters are 
scintillator based, whereas the plug and forward calorimeters are a sandwich of propor- 
tional tube arrays with lead (PEM) or steel (PHA) absorber, and they are all segmented 
into towers which point back to the nominal interaction point. 

The CEM [11] provides electron and photon energy measurements in the region \rj\ < 
1.1 with resolution as/E = 13.5%/VE sin 9 © 1.5%, where E is measured in units of GeV 
and © indicates sum in quadrature. The CEM is physically separated into two halves, 
one covering rj > and one covering 77 < 0. Both halves are divided in azimuth into 24 
wedges that subtend 15° each. Each wedge extends along the z-axis for 246 cm and is 
divided into ten projective towers of approximately 0.1 units in rj. The CEM is 18 radiation 
lengths thick and consists of 31 layers of plastic scintillator interleaved with 30 layers of 
lead sheets. A proportional chamber (CES) measures the electron shower position in the (f) 
and z directions at a depth of ~ 6 radiation lengths in the CEM. The CES module in each 
wedge is a multi-wire proportional chamber with 64 anode wires oriented parallel to the 
beam axis. The cathodes are segmented into 128 strips perpendicular to anode wires. An 
electron and photon shower typically spans several CES channels in each dimension. When 
CTC tracks made by electrons from W boson decays are extrapolated to the CES (r 
184 cm), the CTC extrapolation and the CES shower position match to 0.22 cm (rms) in 
azimuth and 0.46 cm (rms) in z. Both CES/CTC position matching and the CES shower 
shape are used as electron identification variables. 
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The PEM provides energy measurement in the range 1.1< <2.4 and the FEM covers 
2.2< |?7| <4.2. The towers subtend approximately 0.1 in pseudorapidity by 5° in cf). Details 
of the plug and forward calorimeters can be found in [12, 13]. 

All the calorimeters are used to measure missing transverse energy and the central elec- 
tromagnetic calorimeter provides the energy measurement for the electrons in this analysis. 

C: Muon Systems 

Three systems of scintillators and proportional chambers are used to identify muons in this 
analysis. A four-layer array of drift chambers, embedded in each wedge directly outside of 
the CHA, form the central muon detection system (CMU) [14, 15]. The CMU covers the 
region \'q\ < 0.6 and measures a four-point trajectory (called a "stub") with an accuracy 
of 250 /xm per point in r — ^. Charge division gives an accuracy of 1.2 mm per point in 
z. A 0.6 m-thick layer of steel separates the CMU from a second four-layer array of drift 
chambers (CMP). Requiring a muon to have a stub in the CMP reduces the background 
due to hadrons and in-flight decays by approximately a factor of ten. The CMU covers 
approximately 84% of the solid angle for \rj\ < 0.6, while 63% is covered by the CMP, and 
53% by both. Additional four-layer muon chambers (CMX) with partial (70%) azimuthal 
coverage lie within 0.6 < |?7| < 1.0. 

D: Trigger Requirements 

The CDF trigger [16] is a three-level system that selects events for recording to magnetic 
tape. The first two levels of the trigger consist of dedicated electronics. At Level 1, electrons 
arc selected by requiring the presence of deposited energy above 8 GeV in a trigger tower 
(one trigger tower is two physical towers, with a width in pseudorapidity of Ar/=0.2). Muons 
are selected by requiring the presence of a track-stub in the CMU or CMX and, where there 
is coverage, a track-stub in the CMP in coincidence with the CMU. The Level 2 trigger starts 
after a Level 1 trigger has accepted an event. Trigger towers in the calorimeters are combined 
into clusters of total or electromagnetic energy by a hardware cluster finder. Clusters 
and stubs are then matched to tracks found in the CTC by the fast hardware tracking 
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processor. The third-level trigger uses software based on optimized offline reconstruction 
code to analyze the whole event. 
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FIG. 3: One quarter of the CDF detector. The detector is symmetric about the interaction point. 
This is the configuration for Run lb. 

III. DATA SELECTION 

The data presented here were collected by the CDF detector at the Tevatron collider between 
1994 and 1995 ("Run lb"). The signature for a W ^ i + ly event is a lepton with high 
transverse momentum and large missing transverse momentum in the event, due to the 
undetected neutrino. In the electron channel, we select candidate events with the primary 
lepton in the CEM. In the muon channel, the lepton candidate is required to have stubs in 
the CMU, CMP or CMX. These conditions specify what is referred to here as the "central 
lepton" sample. Two samples oi Z ^ e'^ + e~ and Z ^ /j,'^ + are also used for tuning 
the simulation. The details of the trigger requirements can be found in [17]. The integrated 
luminosity is 80±4 pb"-*^. 
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The missing transverse momentum is inferred from the energy imbalance in the event. 
For this purpose, a recoil-energy vector u is defined as the vector sum of the transverse 
energies of all calorimeter towers (including both electromagnetic and hadronic, up to \r]\ < 
3.6), except the ones identified as part of the electromagnetic clusters associated with the 
primary leptons: 

u= ^ EiSmOihi, (4) 

i not i 

where fii is a transverse unit vector pointing to the center of each tower and sin^j is 
computed using the 2;-vertex closest to the electron track, or using the electron track zq 
if there is no z-vertex within 5 cm of the electron track. The vector is a measure of 
the calorimeter's response to jets and particles recoiling against the W. Thus, the missing 
transverse energy (identified with the transverse momentum of the neutrino) is derived as 
f^T = '^{Pt + '^)i where ij- denotes the muon transverse momentum {pt) or the electron 
transverse energy (Et)- The modulus (n) of the recoil vector is an estimator of the W 
boson transverse momentum and it is used to select different ranges of the W boost. 

The analysis uses the transverse mass (My), which is analogous to the invariant mass 
except only the transverse components of the four-momenta are used. Mt is determined 
from the data as 

Mt = ^2P^^t(1-cosA(/)^'^), (5) 

where A^^'' is the angle in the r — (p plane between the transverse momentum of the lepton 
and the missing energy. 

Several selection criteria are chosen to isolate a sample of well measured electrons and 
muons and reduce the backgrounds. For candidates in the W ^ e + u sample, we select 
electrons with Et > 25 GeV and with thepT of the associated track greater than 15 GeV/c. 
Events are accepted only if ^t > 25 GeV. We require a well measured track (crossing all 
eight super-layers of the CTC and with more than 12 stereo hits attached). To exploit the 
projective geometry of the CDF detector, the event vertex reconstructed with the VTX 
is selected to be within 60 cm in z from the origin of the detector coordinates. Fiducial 
requirements are applied to ensure that candidates are selected in regions of well understood 
efficiency and performance of the detector. To remove Z-boson events from the W sample 
a search is made for a partner electron in the central (GEM), plug (PEM), or forward 
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(FEM) calorimeter. Partner electrons are sought with cluster transverse energies greater 
than 20 GeV, 15 GeV and 10 GeV in the GEM, PEM and FEM respectively. Tracks with 
transverse momentum greater than 10 GeV/c and opposite sign to the primary electron are 
also considered. The event is rejected if the invariant mass of the primary electron with 
the partner electron exceeds 60 GeV/c^. The event is also rejected if the partner electron 
is pointing to any non-fiducial volume of the calorimeter, as this may cause the cluster's 
energy to be mis-measured and consequently cause the invariant mass rejection to fail. 

In order to improve electron identification, additional variables are used. They are the 
ratio of the hadronic to electromagnetic deposited energies (Ehad/Eem < 0.1), the match 
between the extrapolated track and the measured position at the GES {Azces < 5 cm), the 
transverse GES shower shape [18] (x^ < 10), and the track isolation (/50o.25 < 1 GeV/c). 
The track isolation variable ISOr is defined as the total transverse momentum from tracks 
(unconstrained by the vertex position) of > 1 GeV/c, that lie within a cone of semi 
opening R = -^(Ar/)^ -|- (A^)^ centered on the lepton track and within 5 cm of the lepton 
z vertex. 

For candidates in the W ^ + u sample, the muon px and the I^t in the event are 
required to be greater than 25 GeV. The quality requirements on the tracks are the same as 
for the electrons. In addition, there are requirements on the impact parameter of the track 
(|cZo| < 0.2 cm) and on the opening angle (> 10°) with the second high-p^ track to remove 
cosmic rays. The muon identification is based on the presence of track-stubs in the muon 
systems and on the deposited energy of the candidates in the calorimeters. The deposited 
energy is required to be less than 2 GeV in the GEM and 6 GeV in the GHA. Furthermore, 
we require that the GTG track, extrapolated at the center of the muon chambers, and the 
track-stub reconstructed in the muon systems match to within 2 cm in the GMU or 5 cm 
in the CMP and GMX. The track isolation cut has not been applied to muon candidates 
since the muon sample is smaller in size and we have preferred a looser selection. The Z 
removal rejects events where there is a second highest-pT (> 10 GeV/c) track in the GTG, 
of opposite sign to the jjl candidate and back-to-back in space (within 10°), that has an 
invariant mass with the /x candidate greater than 50 GeV/c^. 

The Z samples are selected with the same W selection criteria, except the is replaced 
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with a partner high-pT lepton, and the Z removal requirements are not apphcd. Moreover, 
the sample of Z — > e"*" + e~ used for the tuning of the simulation has two CEM electrons, 
both passing electron ID cuts. This choice removes almost all of the QCD background. 

A summary of the selection requirements and the number of surviving events is shown 
in Tables 1 (electrons) and 2 (muons). The accepted samples consists of 22,235 W ^ ^ + 1/ 
candidates and 41,730 W ^ e + u candidates, divided in four recoil ranges. 



Criterion 


W PVPTits Rftpr TPmiirPTTiPTits 


x±x± u J-CAjx o cox l^x\^ 


105 073 


Fiducial requirements 


75,135 


Good electron track 


68,337 


> 25 GeV 


64,254 


El^ > 25 GeV 


54,409 


u < 100 GeV 


54,300 


>15 GeV/c 


52,573 


Mt = 50- 100 GeV/c2 


51,077 


Electron ID 


42,882 


Z removal 


41,730 


u<10 GeV 


31,363 


10 < u < 20 GeV 


7,739 


20 < u < 35 GeV 


2,033 


35 < M < 100 GeV 


595 



TABLE 1: Set of requirements applied to select the W ^ e-\- v data sample. 



IV. MEASUREMENT METHOD 

Ideally one would like to fit the distribution of cos Ocs for the coefficients ai and a2 of 
Eq. (1). However, since the neutrino coming from the W decay is undetected, the kinematics 
of the decay are not completely reconstructed and it is not possible to perform a boost 
into the W rest frame and uniquely determine cos 6'cs- The finite width of the W boson 
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Criterion 


W events after requirements 


Initial ample 


60,607 


E^^^ < 2 GeV and E^^^ < 6 GeV 


56,489 


Not a cosmic candidate 


42,296 


Impact parameter do < 0.2 cm 


37,310 


Track-muon stub match 


36,596 


Good muon track 


33,887 


PT > 25 GeV/c 


29,146 


E!^ > 25 GeV 


25,575 


u<70 GeV 


25,493 


Z removal 


22,877 


Mt=50-100 GeV/c^ 


22,235 


u < 7.5 GeV 


13,813 


7.5 <« < 15 GeV 


5,910 


15 < u < 30 GeV 


2,088 


30 < 7i < 70 GeV 


424 



TABLE 2: Set of requirements applied to select the W ^ ^ + v data sample. 



makes it difficult to solve the equations of the W two-body decay. Even if the mass of 
the W were known on an event by event basis and the detector had perfect resolution, the 
unknown longitudinal component of the neutrino momentum would leave a sign ambiguity 
in determining cos 9c s- 

This measurement therefore exploits the relationship between the transverse mass of 
the W and the lepton polar angle on a statistical basis, i.e. by using the shape of the Mt 
distribution. A similar technique has been successfully applied in [4] to measure 02 from 
W ^ e + u decays. Fig. 4 shows an example of how the distribution of the transverse 
mass of the W changes with different values of a2. Also, since Mt does not contain any 
information on the longitudinal boost of the W boson, it is affected by ai (the forward- 
backward lepton decay asymmetry term) only through residual effects of the geometrical 
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acceptance of the detector. 

The parameter 02 is determined by fitting the Mt distribution to a set of Monte Carlo 
generated templates, each with a different value of a2. A binned log-likelihood method 
is applied to find the best estimate for 0:2 • The procedure is repeated selecting different 
regions of the transverse momentum of the W boson. 




[GeV/c^] 

FIG. 4: Example of the sensitivity of the Mt distribution to a2- Here a2 has been set to and 1, 
and is less than 20 GeV/c. 

V. MONTE CARLO SIMULATION OF W PRODUCTION AND DECAY 

A fast Monte Carlo (MC) generator and a parametrization of the detector response have 
been used in this analysis to simulate W events at CDF [17]. The event generator is based 
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on a leading order calculation of W production and leptonic decay in quark-antiquark anni- 
hilation, including final state QED radiation [19, 20, 21, 22]. The distribution of momenta 
of the quarks is based on the MRS-R2 [23] set of parton distribution functions (PDFs). 
The generated W boson is then Lorentz-boosted, in the center-of-mass frame of the quark- 
antiquark pair, to a specific transverse momentum ■ This measurement uses a broad 
range of pj^, including events at low pj^, where theoretical calculations are not reliable. 
The spectrum of pj^ as a function of the W boson rapidity is therefore derived from p^ 
(the pT oi a Z boson — determined experimentally from Z e"*" + e~, + fi~ events) 
after correcting it by the theoretical prediction for Pj^ /Pt- There is no physics simulation 
of the recoiling jets, instead we model directly the detector response to the recoil against a 
W boson. The parametrization of the detector response and the modeling of the W boson 
recoil up to 20 GeV/c is described in detail in [17]. We have tuned the parameters of the 
model to describe the range of p^ up to 100 GeV/c. Overall, the MC tuning performed for 
this analysis involves: 

— the effects of QCD on the lepton angular distribution, 

— the parametrization of the Z transverse momentum spectrum, up to 
p| = 100 GeV/c, and 

— the detector response to the recoil against high-p^^ Z and W bosons. 



A: Effects of QCD on the lepton angular distribution 

The QCD effects on the lepton angular distribution arc implemented with an event weighting 
procedure in the simulation. Leptons from W decays, generated with a tree-level quark- 
antiquark annihilation, have a purely V — A angular spectrum with a very small distortion 
due to the final state photon emission. Therefore, events are first unweighted by 1/(1 — 
Q cos 9y, where 6 is the lepton polar angle in the parton frame and Q is the lepton charge. 
This effectively factors out any small distortion of the spectrum with respect to a parabola. 
Events are then assigned the appropriate weight (wqcd), where wqcd is defined as a 
function of the lepton angles {Ocs, (t>cs) in the Collins-Soper W boson rest frame: 

wqcd{Ocs, <t>cs) = 1 + cos^ dcs + ^^o(l - 3cos^ ecs) + 
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+- ^2 sin^ 9c s cos 24)cs + 

sin 9c s cos (f)cs - QM cos 9c s- (6) 

Eq. (6) describes the angular modulation induced by the effects of QCD as expressed 
also in Eq. (2), except for the terms with ^1,5,6,7; here they are set to zero, corresponding to 
the Standard Model prediction in the accessible range. The coefficients A2 and Az are 
kept in the angular distribution and assigned the SM dependence with p^, calculated in [3]. 
Notice that the angular coefficients to A2 and ^3 cancel out when integrating analytically 
over 4)cs between and 27r. Nevertheless, detector acceptance effects introduce a small 
residual dependence in the polar-angle spectrum. 

In Eq. (6), wqcd can take negative values if ^0 ^"^^ A4 (or, equivalently, 02 and ai) 
are varied independently in the procedure of fitting for the best parameters. Fig. 5 shows 
the allowed parameter space for a2 and ai. The diagonals in the plot correspond to the 
requirement: 

{1 + a2Cos^ 9cs± Oil cos 9cs) >0, (7) 

for cosOcs = il- The point {ai,a2) = (2, 1) is the Quark Parton Model (QPM) limit in 
the case that the sea quark contribution is neglected, and it has a vanishing cross section 
at &CS = ±180°, as described by the V — A lepton angular distribution. The dotted line is 
the relationship between a2 and ai (at different up to 100 GeV/c), expected from the 
SM including QCD corrections. To prevent wqcd from taking negative values, ai and a2 
are varied only within the allowed region. Note that the sea quark contribution to ai is 
correctly taken into account in the Monte Carlo. 

Because this is an event-weighting procedure, it does not correspond to the inclusion of 
QCD corrections to the generated events: the large-p^^ W events still have to be introduced 
by hand, by imposing a transverse momentum distribution. 



B: Z transverse momentum spectrum 

Prior to the determination of the Z transverse momentum distribution, the Monte Carlo 
simulation is tuned and checked against the Z ^ e'^ + e~ and Z ^ /j,'^ + /j," invariant 
mass distributions from the data. In the electron channel, the Monte Carlo simulation 
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FIG. 5: The ai vs. a2 parameter space. The regions marked with "not aUowcd" are where the 
combination of a2 and ai gives unphysical negative weights to the differential cross section. The 
dotted hne shows the values of ai and a2 at different between and 100 GeV/c. 



reproduces the data with an input Z-mass equal to the world-average [24] within a scale 
factor of 1.0002ib0.0009. In the parametrization of the energy resolution of the CEM: 

aE 13.5% ^ 

K, (8) 



E 

we use K = (1.23 zt 0.26)%. The k term accounts for residual gain variations not corrected 
by the calibration procedure and is obtained from a fit to the Z invariant mass peak. 

There is a small non-linearity correction to extrapolate the energy-scale corrections from 
electrons at the Z pole to the energies typical of aW decay. The average Et for electrons 
coming from Z decay is about 4.5 GeV higher than the Et for W decay. The non-linearity 
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over a small range of energies can be expressed with a slope as 

SEiW) = SE{Z)-[l + ^AET], (9) 

where Se{Z) is the measured scale at the Z pole, ( is the non-linearity factor, and AEt is 
the difference in the average Et between Z and W electrons. The estimate of ^ is derived 
by looking at E/p distributions from the W data and comparing them to the Monte Carlo 
simulation in separate regions of Et- We estimate ^ to be 

^ = -0.00027 ± 0.00005(stat) GeV-^ (10) 

For muons, we use a momentum resolution of 

a{l/pT) = (0.097 ± 0.005) x 10"^ (GeV/c)-\ (11) 

and the reconstructed Z mass peaks in the data and MC match with a ratio of central values 
of l.OOOSibO.OOll. With these inputs, the Monte Carlo simulation reproduces correctly the 
peak position and width of the invariant mass distribution of electron and muon pairs from 
Z bosons. 

Since the QCD corrections to Z production are not included in the Monte Carlo sim- 
ulation, the transverse momentum of the Z bosons needs to be determined from data. 
The distribution is generated in the Monte Carlo simulation using the following ad hoc 
four-parameter functional form: 

fipr) = f^^^ [(1 - Pi)Pi'^'e-^^^ + P.Pi'^'e-^-^^] ; x = p|/(50.0 GeV/c). (12) 

The parameters -Pi,.., 4 are determined from a fit to the observed distribution and 
then corrected to account for the difference between the observed and the generated p^ 
spectrum. Since the difference between the two spectra is very small, the unfolding of the 
effect of the reconstruction is obtained by considering the ratio between them, as predicted 
by the detector simulation. We determine the p^ distribution using separately Z ^ ^'^+fi'' 
and Z ^ e'^ + e" data, and the average is used as the spectrum that is input to the 
Monte Carlo simulation. The uncertainty on the average is used to evaluate the systematic 
uncertainty due to the transverse momentum spectrum determination. Fig. 6 shows the p^ 
distribution of Z ^ n'^ + /j." and Z ^ e'^ + e~ data. The spectra are compared with 
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the simulation where the parameters have been fit to the data. There is a good agreement 
between data and Monte Carlo simulation and the values, normalized per degree of free- 
dom, are very close to 1. 
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FIG. 6: Distributions of from Z — > /i+ + data (a) and Z — > e+ + e data (b) compared with 
the simulation. 



C: Detector response to the recoil against high-pr Z and W bosons 

An estimate of the W boost in the transverse plane comes from the measurement of the 
calorimeter response to jets and particles recoiling against the W. The definition of the 
recoil-energy vector u is given in Eq. (4). The modeling of u in terms of the W boson 
transverse momentum is called the "recoil model" and it is implemented in the Monte 
Carlo simulation of the event. The recoil model is derived using the observed recoil against 
Z bosons, whose kinematics are completely determined by the two leptons. The assumption 
is made that the recoil against Z bosons can be extended to model W events, since the W 
and Z bosons share a common production mechanism and are close in mass. We summarize 
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below the key elements of the recoil model and show how the simulation describes the data 
after fitting the model's parameters to the high-pr ■Z^ boson data. 

1. Recoil Model 

The direction of measured from the reconstructed decay leptons and the perpendicular 
to it form a base in the r — (f) plane on which the recoil vector u can be projected: u= («||, 
u±). The values of u\\ and u± are functions of (addressed here as "response functions") 
with a certain smearing. The smearings are to a good approximation Gaussian distributions 
[5], so that and u± can be parametrized as Gaussians with variable mean and width: 

2. Response functions 

The response function /y is well described by a second order polynomial in the Z transverse 
momentum measured from the reconstruction of the decay leptons. The parameters for 
) are obtained from a fit to — e+ + e"~ and Z — /x+ + /i~ data and the function 
is corrected for a small difference between the true p§, and the observed p§, — which is 
measured from the two leptons' momentum vectors — to feed the correct parameters to 
the simulation. Fig. 7(a) shows the average of which is the response function for the 
parallel component, together with the simulation after fitting for the parameters of f\\. 

is on average smaller than pj^, due to the gaps in the calorimeter and inefficiency in 
the reconstruction of the total energy deposited. Nonetheless, measuring provides an 
estimate of p^ (or ultimately Pj^)- 

The response function f±{pT) is consistent with zero within the statistical uncertainty, 
as expected since u± is the recoil projection perpendicular to p^. The average of u± is 
shown in Fig. 7(b). 

3. Resolutions 

The resolution of the recoil vector components depends on the underlying event and the jet 



18 



activity, au and a± are parametrized in the form: 



omUY^ET) X I |, (14) 



where P2,|| ^^^d P2,± are second order polynomials in p^, whereas cr^b^ contains the un- 
derlying event contribution and is modeled by minimum bias events. In Eq. (14), Umhs is 
expressed as a function of the total transverse energy J^^t, defined as the scalar sum of 
tower transverse energies: 

Y^Et= ^ E^sinOi. (15) 

i not ^± 

Et is a measure of the total transverse energy in the event from all sources, exclud- 
ing the primary lepton. The functional dependence of (Tmbs versus Y^Et is calculated in 
[17]. The explicit dependence in the polynomials is derived here from Z data, using 
both electrons and muons. The parameters are then corrected for the dependence of the 

observed versus the true pj,, as done for the response functions. Fig. 7(c), (d) show the 
resolution of u\\ and u_\_. The resolution a{u\\) worsens at higher p^^ due to increased jet 
activity in the event. The agreements between data and Monte Carlo simulation are good 
in all the plots and the x^'s normalized per degree of freedom are close to 1. 



D: W transverse momentum distribution 

To turn the pj- distribution into a p^ distribution, the simulation applies two weighting 
functions. The first allows for the fact that the pf distribution (as in Eq. (12)) is derived 
with a fit performed to data averaged over all rapidity values (with mean |y|=0.3). However, 
W events need to be generated differentially in both px and y. This weighting function is 
taken from a theoretical calculation of -^^/ (h^) [!'''] • 

The second weighting function turns the pj. distribution, generated with both pr and 
y dependence, into a distribution for the transverse momentum of the W boson. This is 
obtained from the theoretical calculation of g^g^ ^/ ^ [25, 26, 27, 28]. Resummed 

calculations are used for correcting the difference between the W and the Z pr distributions. 
The ratio is between 0.9 and 1.0 over the pr range of interest. Since this is a ratio, the 
uncertainty is expected to be small because of cancellation of systematics. Indeed, by 
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FIG. 7: (a) and (b) Comparison of the data with the simulation for the recoil response components 
Mil and u± versus . (c) and (d) the resolutions cr{uu) and a{u±) versus . 
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varying the PDF, a^, or the type of calculation, the resulting uncertainty in pj. is small in 
comparison to the uncertainty arising from the statistics of the Z sample used to define the 
distribution [29, 30, 31, 32]. 

Although due to the undetected neutrino we cannot compare directly the spectrum 
in the simulation with the data. Fig. 8 shows a comparison of the recoil against the W 
in the electron and muon channel. The recoil includes the distribution as well as all 
the response and resolution parameters derived using the Z sample. The shaded band 
corresponds to the uncertainty on the spectrum only. Since the recoil model and the 
Pj^ spectrum are derived with a sample that is much smaller than the W sample, there is 
a degree of freedom in optimizing the parameters to improve the agreement with W data. 
However, we choose not to optimize the parameters using any of the W boson distributions 
to prevent a possible source of bias when fitting the transverse mass distribution. We treat 
the statistical uncertainty of the recoil model and pj^ spectrum as a source of systematic 
uncertainty for a2- 




FIG. 8: Distribution of the recoil against the W boson compared with the simulation 'mW~^e + i^ 
data (a) and W ^ ^jl + u data (b). 
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VI. BACKGROUNDS 



There are three main sources of background to the W ^ i + v data sample of this analysis 
(where £ stands either for an electron or a muon): 

• W— > r + 1/ events, with the r subsequently decaying into a muon or electron and two 
neutrinos. 

• Z— >■ + i~ events, where one of the leptons is not detected. 

• QCD dijet events, where a jet is wrongly identified as a lepton and the total energy 
in the event is incorrectly measured to give a signal. 

There is a small background contribution from tt decays, which is estimated to be ~25 
events in the electron channel and ~12 in the muon channel [33] and affects the high recoil 
range only. The background from cosmic rays in the muon channel is approximately 0.2% 
[17] of the total W ^ ^ + u candidates, with a flat Mt distribution. This corresponds to a 
negligible contribution compared with the dominant backgrounds. 

A shape for the transverse mass distribution is determined for each background source 
and added to the transverse mass distribution of the simulated W events. For background 
the shape is taken from [34]. 

A: W ^ T -\- u background 

The background from W ^ t + v events, where the r decays leptonically, is virtually 
indistinguishable from the VF— > e + v ot W ^ n + v signal. The event generator used 
for the simulation of W events in this analysis is capable of simulating W ^ t + u, where 
the r lepton is then decayed into /x + 2i/ or e + 2i/. The background level is found to be 
approximately 2% of the total W sample, with softer charged lepton pT and spectra. 
The W ^ T + v background fractions are listed in Tables 3 and 4 for the electron and muon 
channel respectively. The shape of the transverse mass distribution is also taken from the 
Monte Carlo simulation oi W ^ t + v events, separately for each of the W boson recoil 
ranges. 

B: Z ^ e+ +£- background 
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Z events enter the W sample when one of the leptons is not detected ( "lost leg" ) and there 
is missing transverse energy in the event. 

1. Electron channel 

As part of the W candidate selection procedure the primary electron is always required to 
have been detected in the central calorimeter. The Z removal procedure ensures the rejec- 
tion of events with a second oppositely charged high-p^ track, or high-energy calorimeter 
cluster, and invariant mass of the electron-candidate pair compatible with a Z boson decay 
(Mee > 60 GeV/c^). When the track associated with the second electromagnetic cluster is 
pointing to any non-fiducial volume of the calorimeter, the event is rejected irrespective of 
the invariant mass value. This ensures that the event would still be rejected if the second 
electron has emitted a photon and the invariant mass with the primary electron track falls 
outside the Z invariant mass exclusion range. 

The Monte Carlo simulation is used to estimate the Z background due to the inefficiency 
of the calorimeters in detecting the second leg, or when the second electron points beyond 
the coverage of the forward calorimeter (jr/l > 4.2). The total background level from Z 
events in the electron channel is very small, and is listed in Table 3. 









Recoil [GeV] 






Type: 


(0-10) 


(10-20) 


(20-35) 


(35-100) 


All 


W ^T^V 


2.15 


1.74 


1.31 


1.57 


2.01 


Z ^e + {e) 


0.00 


0.02 


0.12 


0.39 


0.01 


QCD jets 


0.23±0.11 


0.39±0.14 


0.14±0.10 


0.5±0.3 


0.26±0.12 


ti 


0.00 


0.00 


0.49±0.20 


2.50zb0.80 


0.06±0.02 


Total 


2.38±0.11 


2.15±0.14 


2.06±0.22 


4.96±0.85 


2.42±0.12 



TABLE 3: Summary of the backgrounds to — > e -|- (as percentages of the W candidate sample) 
in different W recoil ranges. The uncertainty is negligible fov W ^ t + u and Z — > e -|- (e). 

2. Muon channel 
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The event seleetion applied in this analysis removes events with opposite sign tracks (found 
in the CTC) that combine with the identified muon to give an invariant mass greater than 
50 GeV/c^. The number oi Z fi'^ + fi" events not removed by the Z selection criteria is 
consistent with zero when both muons pass through the fiducial tracking volume {\rj\ < 1). 

However, a significant number of Z events may enter the W sample when one of the 
muons goes outside the fiducial tracking volume. About 20% of Z ^ /j,'^ + n" events have 
one of the muons outside \r]\ < 1, either at the edge of the tracking volume (|ry| ~ 1.1) or at 
higher rj, beyond the coverage of the CTC. The estimate of the background in these cases is 
based on the simulation, which uses the tracking efficiency map determined using electrons 
detected in the calorimeter from the W ^ e + v sample. The background level found is of 
the order of 4% and it is listed in Table 4. The shape of the transverse mass distribution 
of lost-leg events is also derived from the Monte Carlo simulation. 



Type: 


(0-7.5) 


(7.5-15) 


Recoil [GeV] 
(15-30) 


(30-70) 


All 




2.24 


1.94 


1.63 


2.37 


2.11 


Z ^ + 


4.25 


4.00 


3.67 


2.95 


4.11 


QCD jets 


0.45±0.19 


0.79±0.29 


0.81±0.52 


1.40±1.18 


0.59ib0.26 


tt 


0.00 


0.00 


0.19±0.09 


1.89±0.70 


0.05±0.02 


Total 


G.94±0.19 


6.73±0.29 


G.30±0.53 


8.61±1.37 


G.86±0.26 



TABLE 4: Summary of the backgrounds toW ^ n + f {as percentages of the W candidate sample) 
in different W recoil ranges. The uncertainty is negligible fov W ^ t + u and Z — > + {/j.). 

C: QCD background 

Dijet events can pass the W selection cuts if one of the jets is mis-identified as a lepton and 
one of them is incorrectly measured and gives a high missing- E't signal. This is referred to 
as QCD background. W candidate events which arc background from QCD would typically 
have the charged lepton or the neutrino predominantly back-to-back or collinear with the 
leading jet. Real W events, on the other hand, have a nearly uniform distribution of the 
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lepton-jet opening angle, at least for low . For higher , W events also exhibit a slight 
tendency to have the leading jet, which is recoiling against the W , in the opposite direction 
to the charged lepton and the neutrino. 

1. Electron channel 

Fig. 9(a) shows the distribution of the opening angle in the r — 4> plane between the electron 
and the leading jet. The leading jet is the highest energy jet in the event with energy of 
at least 5 GcV. The plot shows three samples enriched in QCD background together with 
the W candidates sample. Two of the enriched QCD samples are derived by reversing the 
electron ID cuts in the W preselection sample. The third is taken from dilepton events {Z 
candidates that do not pass the opposite charge requirement on the two leptons) which we 
refer to as the QCD control sample. The samples enriched in QCD all show the expected 
peaks at 0° and 180°. 
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FIG. 9: (a) Azimuthal angle between the electron candidate and the leading jet in the QCD samples 
and in the W-candidates sample, (b) Number of events in the plane of recoil versus isolation in the 
QCD-enriched sample, derived from the dilepton sample with a same-sign requirement and all the 
electron-ID cuts applied. 
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When the W recoil is less than 20 GeV the background is estimated by counting the 
excess of events in the distribution of A(p{i — jet). The signal component is estimated by 
fitting a linear function to the middle part of the A<p{£ — jet) distribution. Almost all the 
W candidates with recoil greater than 10 GeV come associated with at least one jet, and 
we account separately for events that do not have an associated 5 GeV jet. Since the 
W candidates greatly outnumber the background events when the electron is isolated, the 
counting is done in bins of increasing isolation, and the background is extrapolated back to 
the signal region of ISOq,25 < 1 GeV/c. The same background estimate is cross-checked 
by selecting events at high isolation (6 < ISO0.25 < 10 GeV/c) and using the fraction of 
isolated to non-isolated QCD events, seen in the QCD control sample, to predict the number 
in the signal region. Fig. 9(b) shows the two-dimensional distribution of the recoil versus 
lepton isolation in the QCD control sample. 

We estimate 74ib36 background events due to QCD in the 0—10 GeV recoil range and 
SOibll in the 10—20 GeV recoil range. This includes an additional 10±7 events in the 
0—10 GeV recoil bin due to W events with no leading 5 GeV jet, as derived from the 
fraction of events with and without a jet in the QCD control sample. The uncertainties 
include a systematic component due to the method. At higher W recoil the estimate of 
the background is 3=b2 events in both the 20—35 and 35—100 GeV bins. This is estimated 
with both the QCD control sample (by using the ratio of low to high recoil) and the direct 
counting of the excess of events at 0° and 180°. In the latter, the non-uniform opening 
angle distribution of the recoiling jet and W^-decay leptons is partially accounted for by a 
slope in the fit to the opening angle distribution. The small background contribution makes 
it unnecessary to accurately model the signal angular distribution. 

The shape of the transverse mass distribution of the QCD background is obtained by 
reversing the isolation cut and selecting events with anti-isolated electron tracks. The Mt 
distribution shapes, at different recoil ranges, are seen to be largely independent of the 
anti-isolation cut. Fig. 10 shows the Mt distribution of the backgrounds in the electron 
channel, scaled by the estimated amount as a percentage of the W candidates. 

2: Muon channel 

QCD events can mimic W— fi+u mainly in two ways. The first is when a heavy flavor quark 
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in one of the jets decays into particles that include a high-pT muon (e.g. b ^ c + iJ, + v). 
In order for the mnon and neutrino to have enough pT to pass the W selection cuts, the 
6 quark needs to have a high transverse momentum, which leads to small opening angles. 
Therefore this type of event will have the muon and the neutrino almost parallel to one 
of the jets. The second major type of QCD background process occurs when a hadron is 
misidentified as a muon. The energy of one of the jets should also be incorrectly measured, 
in order to give the appearance of a high missing- i^y signal. In this case, the neutrino and 
the muon will be reconstructed either nearly parallel to one jet or back-to-back and parallel 
to the two jets. Moreover, in both the processes considered, the muon is not likely to be 
isolated. 

The QCD background to W ^ fi + u events is estimated in the same way as for the 
electron channel in the four recoil bins. We expect 62 it 26, 47ibl7, 17ibll, and 6ib5 events 
in the four recoil ranges. Fig. 11 shows the Mt distribution of the backgrounds in the muon 
channel scaled by the estimated amount as a fraction (percent) of the W candidates. 



A set of Monte Carlo generated templates of the Mt distribution is compared to the dis- 
tribution derived from the data. When each template distribution is compared to the data, 
a likelihood is computed according to: 



where the sum runs over the number of bins of the Mt histogram, nf is the number of 
entries in each bin of the data histogram, and pf^^ are the probabilities per-bin. The values 
of ^^'-'(0:2) = «i^*"/"'toF given by the entries in the template histogram, one template 
for each value of 0:2 ■ The maximum of the likelihood function locates the best estimate 
for the value of 0:2 ■ Fig- 12 shows the likelihood functions in four different regions for 
the electron and muon channels. The likelihood functions have been shifted vertically so 
that the maximum is always at zero. The la statistical uncertainty on each fit is evaluated 



VII. FITS AND SYSTEMATIC UNCERTAINTIES 



A: The likelihood fits 




(16) 



i=l 
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FIG. 10: Electron channel: the transverse mass distribution from the background sources in four 
W recoil ranges. The plots are in percentage of the W data in the specific region. 
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FIG. 11: Muon channel: the transverse mass distribution from the background sources in four W 
recoil ranges. The plots are in percentage of the W data in the specific region. 
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at the points on the hkehhood curve which arc 1/2 unit below the maximum. The four 
recoil regions are 0—10, 10—20, 20—35, and 35—100 GeV/c for the W e + v data and 
0-7.5, 7.5-15, 15-30, and 30-100 GeV/c for \h.e W ^ ji + v data. The choice of the 
ranges is constrained by the sample size in the high-p^ regions, due to the rapidly falling 
distribution. Moreover, the smaller sample of the muon channel is reflected in the recoil 
ranges covering lower values than in the electron channel. Tables 5 and 6 summarize 
the results of the fits for 02- Fig. 13 and 14 show the transverse mass distribution of the 
data compared with the simulation, where a2 has been set to the best-estimate values. 

B: Systematic uncertainties 

Systematic uncertainties on the measurement of a2 for this analysis derive from the simu- 
lation of W events, the detector response, and the estimate of backgrounds. Some of these, 
although classified as systematic, may be statistical in nature. This is the case for the 
detector recoil response and the W transverse momentum spectrum, since they arc derived 
from the Z ^ e'^ + e'^ and Z ji'^ + jjL'^ data samples. In the following, each source 
of systematic uncertainty is discussed and an estimate is determined for the shift on the 
measured values of a2- Tables 5 and 6 contain a summary of the various contributions and 
the total systematic uncertainty. 

1. Event Selection Bias 

The electron isolation requirement may introduce a bias on the measurement of 02- For 
example, if the electron travels close to the recoil, there is greater opportunity for the event 
to be rejected. Also, there could be a correlation of the selected sample with a2, which is 
correlated with the QCD activity in the event. This bias is investigated by removing the 
isolation requirement, evaluating appropriately the increase in backgrounds, and measuring 
the change in a2- The maximum shifts observed are within the systematic uncertainty of 
the background determination. Moreover, by changing widely 02 in the simulation, the 
spectrum of the opening angle between recoil and electron directions is not significantly 
affected. We do not apply an isolation requirement to the muon channel. 

2. Parton Density Functions 
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FIG. 12: Likelihood functions of the fits for oli, in the four PF-boson recoil regions for the electron 
and muon channels. 



31 



3000 



> 



& 

> 



2000 - 



1500 



1000 - 



500 







. W— > e+v 


• CDF Run lb data 


u: 0-10 GeV 


— Monte Carlo 




X"/dof = 34/25 


-♦4- 









50 55 60 65 70 75 80 85 90 95 100 

2 

Transverse Mass [GeV/c ] 



a 



> 
o 

& 

> 



700 
600 
500 
400 
300 
200 
100 







L e+V 


• CDF Run lb data 


u: 10-20 GeV 


— Monte Carlo 




i X^/dof= 22/25 




^„ , 1 , , , , 1 , , , , 1 , , , , 1 , , , 



50 55 60 65 70 75 80 85 90 95 100 

2 

Transverse Mass [GeV/c ] 



(b) 




FIG. 13: Comparison of the transverse mass distribution from the W ^ e + v data (filled circles) 
with the simulation (solid line) in four recoil regions. In the Monte Carlo simulation, a2 has been 
set to the best-fit value for each recoil range. The shaded histograms indicate the background 
contribution that is estimated to be present in the data and that has been added to the simulation. 
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FIG. 14: Comparison of the transverse mass distribution from the W ^ ji + v data (filled circles) 
with the simulation (solid line) in four recoil regions. In the Monte Carlo simulation, a2 has been 
set to the best-fit value for each recoil range. The shaded histograms indicate the background 
contribution that is estimated to be present in the data and that has been added to the simulation. 
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Recoil range [GeV] 


0-10 


10-20 


20-35 


35-100 


Q!2 measured 


1.09 


1.14 


0.67 


-0.22 


Statistical uncertainty 


±0.05 


±0.13 


±0.29 


±0.36 


a2 predicted 


0.98 


0.84 


0.55 


0.25 


Mean [GeV/c] 


6.2 


15.9 


33.3 


59.2 


Nevt 


31363 


7739 


2033 


595 


Systematic uncertainties: 










PDFs 


±0.01 


±0.01 


±0.01 


±0.01 


W mass 


±0.02 


±0.01 


±0.04 


±0.04 


Input pj. 


±0.02 


±0.03 


±0.03 


±0.04 


Recoil model 


±0.01 


±0.05 


±0.04 


±0.20 


Backgrounds 


±0.01 


±0.01 


±0.01 


±0.01 


C<)ml)iii(Hl syst ic 


±0.03 


±0.0C 


±0.07 


±0.21 



TABLE 5: Summary of the measurement of 0:2 with the W ^ e + u data. The mean pj!^ corre- 
sponding to each recoil range is the mean of the distribution of the "true" W transverse momentum 
in the Monte Carlo simulation. 

The parton distribution functions are used in the Monte Carlo simulation to determine 
the quark content of the proton, and hence the rapidity distribution of the generated W 
bosons. The set of PDFs used to simulate the events in this analysis is MRS-R2 [23]. These 
PDFs describe well the CDF low-ry t^-charge asymmetry data. To evaluate the impact 
of the choice of PDFs on the measurement of 02, two Monte Carlo samples have been 
generated with MRMS-D~ and MRMS-DO, sets that were not tuned on CDF data and 
differ significantly from MRS-R2. 02 has been measured with both sets of PDFs. The 
observed shifts are ±0.01 in all recoil regions, a small fraction of the statistical uncertainty. 
This is conservatively taken to be the systematic uncertainty due to the choice of PDFs. 

3. The W mass value 

The transverse mass distribution is sensitive to the value of the W mass used in the Monte 
Carlo simulation. The dependence comes from the fact that the transverse mass spectrum 
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Recoil range [GeV] 


0-7.5 


7.5-15 


15-30 


30-70 


a2 measured 


1.03 


1.24 


0.74 


0.24 


Statistical uncertainty 


±0.08 


±0.18 


±0.40 


±0.51 


ai predicted 


0.99 


0.92 


0.70 


0.32 


Mean [GeV/c] 


5.4 


11.1 


24.7 


49.7 


Nevt 


13813 


5910 


2088 


424 


Systematic uncertainties: 










PDFs 


±0.01 


±0.01 


±0.01 


±0.01 


W mass 


±0.02 


±0.01 


±0.04 


±0.04 


Input pj. 


±0.02 


±0.03 


±0.03 


±0.04 


Recoil model 


±0.01 


±0.05 


±0.04 


±0.20 


Backgrounds 


±0.01 


±0.02 


±0.02 


±0.03 


C( )ml )in(Hl syst (^ual ic 


±0.03 


±0.0G 


±0.07 


±0.21 



TABLE 6: Summary of the measurement of 0:2 with the W ^ fj, + v data. The mean corre- 
sponding to each recoil range is the mean of the distribution of the "true" W transverse momentum 
in the Monte Carlo simulation. 

has a Jacobian peak at about the value of the W mass. The value of the W mass in 
the Monte Carlo simulation is set to the LEP average [35] 80.412±0.042 MeV/c^, in order 
to be independent of the value measured at CDF. An uncertainty on M\y of 40 MeV/c^ 
corresponds to a systematic uncertainty on a2 of 0.01—0.04. 

4- spectrum 

The W transverse momentum spectrum is derived from the Z sample by measuring p^, 
and using the relatively well known ratio of the Pj^ /p^ distributions from theory. The 
Pj^ distribution is measured from both the Z ^ + e~ and Z /i^ ± fi~ data, and 
then averaged. To account for statistical and systematic uncertainties in determining the 
p^ spectrum, additional MC datasets are generated using the p^ from the electron or the 
muon Z-decay channels only. The measured a2 shifts by between 0.02 and 0.04. 

5. Recoil Model 
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The recoil model consists of response and resolution functions derived from the Z + e~ 

and Z + /x~ data. There are statistical uncertainties in the coefficients of the model, 

which are used here to evaluate a systematic uncertainty. Each of the parameters is changed 
and the 0:2 value is measured. The dispersion of the set of new measurements is taken as 
the systematic uncertainty, which increases with as shown in Tables 5 and 6. The recoil 
model is one of the main sources of uncertainty here since it is constrained with a statistical 
sample much smaller than the W sample itself. The impact of a slight disagreement between 
the W recoil distribution in data and simulation has been estimated by shifting the edges 
of the recoil ranges one at a time by 0.1 GeV/c, only in the data but not in the Monte 
Carlo simulation, simulating event migration between bins. The value of 0.1 GeV/c is the 
difference between the mean of the recoil distributions in the data and in the MC. The 
coefficient a2 has been observed to shift between 0.01 and 0.04 in the four bins. This is 
included in the quoted systematic uncertainty due to the recoil model. 

6. The angular coefficients and ai input value 

Although the distribution of | cos 6 cs\^ a-nd hence M^, should only depend on 0.2 and all the 
remaining angular coefficients should integrate out in practice geometric acceptance causes 
some angular coefficients to have a residual effect on the shape of the Mt distribution. 
Coefficients Ai, A5, Ag, Ay are predicted to be negligible in the Standard Model and are 
set to zero. A2 and A3 are kept in the angular distribution (see Eq. (2)) and are set to 
their Standard Model values. As an estimate of the sensitivity to these terms, neglecting A2 
and A3 results in a shift in the value of 0:2 of 0.02—0.07 in the four bins. These values 
are not included in the systematic uncertainty since the uncertainty on the theoretical SM 
calculation is expected to be much smaller than 100%. 

The coefficient ai also enters the Mt distribution. However, when fitting for 02 at low 
, a\ cannot be set to the SM expected value, due to the requirement of positive event 
weights expressed in Eq. (7). a\ is therefore set to 2y^, which lies in vicinity of the SM 
path for low p^ . With this choice, Eq. (7) translates into a condition for (1 it cos ^05)^, 
which is always positive and prevents assigning negative weights in the region around the 
Quark Parton Model point. A negligible change in the measured 0:2 is visible by setting 
o^iiPr ) to different paths around the SM expectation. For higher p^ (> 20 GeV/c), ai 
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is set to the full SM prediction as there is no danger of assigning negative weights in that 
region. 

7. Backgrounds 

The main sources of uncertainty due to backgrounds come from the estimates of the QCD 
and tt contributions. The QCD background is estimated from the data using the lepton 
isolation and the angular distribution between the lepton and the jets in the event and the 
tt background is taken from [33]. The systematic uncertainty on the measured values of 
a2 is derived by changing the QCD and ti background contents in each pj! range by the 
uncertainty given in the background estimate results in Tables 3 and 4. A maximum shift 
of 0.03 on a2 is observed. 

VIII. RESULTS AND CONCLUSIONS 

Fig. 15 shows the results of this measurement on a plot of a2 versus pjf. The position of 
the points on the x axis has been determined by using the mean of the Monte Carlo distri- 
bution of p}^ corresponding to each recoil range. The solid curve represents the Standard 
Model prediction reported in [3]. The trend is a decrease of 02 with increasing pj^, which 
corresponds to an increase of the longitudinal component of the W polarization. The rate 
measured from a linear fit is ^^15% per 10 GeV/c. The four measurement points from the 
electron channel can be used together with those from the muon channel to compute a 
with respect to the Standard Model expectation. The result is x^=1.5, normalized for 8 
degrees of freedom and considering statistical and systematic uncertainties. 

The measurements of a2 with the electron and muon channels are combined in Fig. 16 
and Table 7. The position in pjf is determined by a weighted mean of each pair of electron 
and muon measurements. The values of 02 are then scaled at the common p}^ value using 
a linear fit and then averaged taking into account the size of the statistical uncertainties. 
Systematic uncertainties are completely correlated between the electron and muon channels. 
The triangles are from [4] and represent the current best values. 

In conclusion, we have measured the leptonic polar-angle distribution coefficient a2 as 
a function of the transverse momentum of the W boson. The results obtained from the 
electron and muon channels are combined together and the measurement reduces by about 
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[GeV/c] 


a2 (CDF combined) 


a2 (theory) 


5.9 


1.07±0.04(stat)±0.03(syst) 


0.98 


13.9 


1.18±0.10(stat)±0.06(syst) 


0.89 


29.7 


0.70±0.23(stat)±0.07(syst) 


0.61 


55.3 


-0.05±0.29(stat)±0.21(syst) 


0.23 



TABLE 7: Summary of the measured a2 combining the electron and muon channels. 

50% the uncertainty on the current best values up to ~ 30 GeV/c. The result is in good 
agreement with the Standard Model expectation up to NLO, whereby 02 decreases with 
as a consequence of QCD corrections to the W polarization. Since the uncertainty is largely 
dominated by statistics especially at higher W transverse momenta, this measurement can 
significantly benefit from the larger data sample pp collisions at -y/s = 1.96 TeV that is 
being collected at CDF in Run II. 
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